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Abstract. An improved code for the calculation of astrophysical reaction 
rates in the statistical model is presented. It includes the possibility to study 
isospin effects. Such effects heavily affect rates involving self-conjugate nuclei 
and may also be found in reactions on other intermediate and heavy targets. 
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^ ■ 1. Introduction 

The investigation of explosive nuclear burning in astrophysical environments is a chal- 
lenge for both theoretical and experimental nuclear physicists. Highly unstable nuclei 
are produced in such processes which again can be targets for subsequent reactions. 
Cross sections and astrophysical reaction rates for a large number of nuclei are required 
to perform complete network calculations which take into account all possible reaction 
links and do not postulate a priori simplifications. 

The majority of reactions can be described in the framework of the statistical 
model (compound nucleus mechanism, Hauser-Feshbach approach, HF) [|TJ, provided 
that the level density of the compound nucleus is sufficiently large in the contributing 
energy window [0] . In astrophysical applications usually different aspects are emphasized 
than in pure nuclear physics investigations. Many of the latter in this long and well 
established field were focused on specific reactions, where all or most "ingredients", 
like optical potentials for particle transmission coefficients, level densities, resonance 
energies and widths of giant resonances to be implemented in predicting El and Ml 7- 
transitions, were deduced from experiments. As long as the statistical model prerequisites 
are met, this will produce highly accurate cross sections. For the majority of nuclei 
in astrophysical applications such information is not available. The real challenge is 
thus not the well-established statistical model, but rather to provide all these necessary 
ingredients in as reliable a way as possible, also for nuclei where none of such information 
is available. 

In Section an improved code for the calculation of astrophysical reaction rates in 
the statistical model is briefly presented. It includes the possibility of studying isospin 
effects. The latter are further discussed in Section |3|. 
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2. The code NON-SMOKER 



Based on the well-known code SMOKER ||, an improved code for the prediction of 
astrophysical cross sections and reaction rates in the statistical model has been developed. 
The current status of the new code NON-SMOKER is described in the following. 

The final quantities entering the expression for the cross section in the statistical 
model Jl| are the averaged transmission coefficients. They do not reflect a resonance 
behavior but rather describe absorption via an imaginary part in the (optical) nucleon- 
nucleus potential [|J. In astrophysics, usually reactions induced by light projectiles (neu- 
trons, protons, a particles) are most important. Global optical potentials are quite well 
defined for neutrons and protons. It was shown |3], ||, || that the best fit of s-wave 
neutron strength functions is obtained with the optical potential by @, based on mi- 
croscopic infinite nuclear matter calculations for a given density, applied with a local 
density approximation. It includes corrections of the imaginary part || ||. A similar 
description is used for protons. 

Optical potentials for a particles are treated in the folding approach flQ| , with a 
parametrized mass- and energy-dependence of the real volume integral [11]. The mass- 



and energy-dependence of the imaginary potential is parametrized according to ]I2| and 
additionally includes microscopic and deformation information [fUjl - 

Deformed nuclei are treated by an effective spherical potential of equal volume |3], ^[ . 
For a detailed description of the formalism used to calculate El and Ml 7-transmission 
coefficients and the inclusion of width fluctuation corrections [13], see and references 
therein. 

The level density treatment has been recently improved . However, the problem 
of the parity distribution at low energies remained. The new code includes a modified 
version of the description @, accounting for non-evenly distributed parities at low ener- 
gies, based on most recent findings within the framework of the shell model Monte Carlo 
method [[U] . 



energies 
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Additionally, the included data set of experimental level information (excitation 
spins, parities) has been updated as well as the experimental nuclear 
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These data bases are continuously updated. For theoretical masses, there is 
a choice between different mass models (e.g. by Hilf et al ]17]], FRDM ETFSI []19fl), of 
which currently the FRDM is favored. Microscopic information needed for the calculation 
of level densities and a+nucleus potentials are also taken from the FRDM, as well as 
experimentally unknown ground state spins |2"0fl . 

Finally, isobaric analog states T > = T < + 1 = T g s - + 1 are explicitly considered in 
the new code. This will be discussed in the next section. 



3. Inclusion of isospin 

The original HF equation |l| implicitly assumes complete isospin mixing but can be 
generalized to explicitly treat the contributions of the dense background states with 



isospin T < = T g f3, and the isobaric analog states with T > = T < + 1 [21], |22|, |23|, |24 
In reality, compound nucleus states do not have unique isospin and for that reason an 
isospin mixing parameter /j, J. was introduced ||21||, which is the fraction of the width of 



T > states leading to T < transitions. For complete isospin mixing /j, I— 1, for pure T < 
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states p 1= 0. In the case of overlapping resonances for each involved isospin, p J. is 
directly related to the level densities p < and p > , respectively. Isolated resonances can 
also be included via their internal spreading width 1^ and a bridging formula was derived 



to cover both regimes [25 



In order to determine the mixing parameter p |= p j (E), experimental infor- 
mation for excitation energies of T > levels is used where available [|T5| , [2~6|j in the code 
NON-SMOKER. Experimental values for spreading widths are also tabulated |26 . 
Similarly to the standard treatment for the T < states, a level density description is 
invoked above the last experimentally known T > level. Since the T > states in a nucleus 
(Z,N) are part of multiplet, they can be approximated by the levels (and level density) 
of the nucleus (Z— 1,N+1), only shifted by a certain energy This displacement en- 
ergy can be calculated J27J and it is dominated by the Coulomb displacement energy 
E d = £Coui + £ _ 

The inclusion of the explicit treatment of isospin has two major effects on statisti- 
cal cross section calculations in astrophysics: the suppression of 7-widths for reactions 
involving self-conjugate nuclei and the suppression of the neutron emission in proton- 
induced reactions. Non-statistical effects, i.e. the appearance of isobaric analog reso- 
nances, are included in the treatment of the mixing parameter p j [gfj but will not be 
further discussed here. 

3.1. 7 -Widths 

The isospin selection rule for El transitions is AT = 0, 1 with transitions 0^0 being 
forbidden [p8 ]. An approximate suppression rule for AT = transitions in self-conjugate 



nuclei can also be derived for Ml transitions [p8 |. 

In the case of (0,7) reactions on targets with N = Z, the cross sections will be 
heavily suppressed because T = 1 states cannot be populated due to isospin conservation. 
A suppression will also be found for capture reactions leading into self-conjugate nuclei, 
although somewhat less pronounced because T = 1 states can be populated according 
to the isospin coupling coefficients. 

In previous reaction rate calculations |29], ||, |3(J the suppression of the 7-widths was 
treated completely phenomenologically by dividing by quite uncertain factors of 5 and 2, 
for (0,7) reactions and nucleon capture reactions, respectively. In the new code NON- 
SMOKER, the appropriate 7-widths are automatically obtained, by explicitly accounting 
for population and decay of T < and T > states, and considering isospin mixing by the 
parameter p [. 

3.2. Competition cusps in proton-induced reactions 

Assuming incomplete isospin mixing, the strength of the neutron channel will be sup- 



pressed in comparison to the proton channel in reactions p+target [23[ P4| . This 
leads to a smaller cross section for (p,n) reactions and an increase in the cross section 
of (p,7) reactions above the neutron threshold, as compared to calculations neglecting 
isospin (i.e. implicitly assuming complete isospin mixing with p {= 1). 

The isospin mixing parameter was varied in the theoretical investigation of a 
51 V(p,7) 52 Cr experiment [BT]]. It was found [RTJ that complete isospin mixing closely 
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reproduced the measured cross sections when width fluctuation corrections were consid- 
ered. Width fluctuation corrections affect the (p,7) cross sections above as well as 
below the neutron threshold, whereas incomplete isospin mixing only reduces the cross 
sections above the threshold. Thus, the two corrections can be discriminated. Mainly 
from that result, it was concluded that — contrary to width fluctuation corrections - 
isospin can be neglected. 

However, a closer investigation of the T> levels in 52 Cr (using p5[ and [H) shows 



that isospin mixing should be rather complete already at the neutron threshold (since 
the first T > state is almost 1 MeV below the threshold This is also true for 

lighter targets. Nevertheless, for reactions on more heavy nuclei (Z > 30) the neutron 
and proton threshold, respectively, will still be in a region of incomplete isospin mixing 
and therefore isospin effects should be detectable there. This effect, however, does not 
play such an important role in the calculation of astrophysical reaction rates as the 
suppression of the 7-width described in the previous chapter, because of the averaging 
over an energy range (the Gamow window) in the calculation of the rate. 



4. Conclusion 

The new code NON-SMOKER makes use of the latest set of descriptions for the calcula- 
tion of the nuclear properties needed to reliably predict astrophysical reaction rates, such 
as masses, level densities, nucleon- and a-potentials, GDR energies and widths, width 
fluctuation corrections. Additionally, the possibility of studying isospin effects has been 
included. This also leads to a more fundamental treatment of the 7-width suppression 
for compound nuclei with T g,s ' = T < = 0. 

Nevertheless, more experimental data are needed to check and further improve 
current parametrizations. Especially investigations over a large mass range would prove 
useful to fill in gaps in the knowledge of the nuclear structure of many isotopes and 
to construct more powerful parameter systematics. Such investigations should include 
neutron-, proton- and a-strength functions, as well as radiative widths, and charged 
particle scattering and reaction cross sections for stable and unstable isotopes. More 
capture data with self-conjugate final nuclei would also be highly desireable. 

This information can be used to make future large-scale statistical model calcula- 
tions even more accurate. 
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